Dynamics of Vortex Nucleation in Nanomagnets with Broken Symmetry 
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We investigate fundamental processes that govern dynamics of vortex nucleation in sub- 100 nm 
mesoscopic magnets. We focus on a structure with broken symmetry - Pacman-like nanomagnet 
shape - in which we study micromagnetic behavior both by means of a simple model and numerically. 
We show that it is possible to establish desired vortex chirality and the polarity of vortex core by 
applying only quasi-static in-plane magnetic field along specific directions. We identify the modes of 
vortex nucleation that are very robust against external magnetic field noise. These vortex nucleation 
modes are common among wide range of sub- 100 nm magnets with broken rotational symmetry. 

PACS numbers: 75.75-c, 75.60.Jk, 75.78.Cd 
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Confinement leads to fundamental changes in physical 
behavior of materials due to increased role of the sur- 
face. In mesoscopic magnetic materials such changes in 
the energy landscape could lead to novel magnetic spin 
configurations such as vortices. Equilibrium properties 
of these new topological states are governed by both the 
properties of the magnetic material and the geometry of 
the object. On the other hand, confinement also leads to 
distinct dynamic behavior of the new topological states, 
since the available energy levels are very much limited. 
The transition probabilities between different states can 
thus be controlled by careful engineering of the geom- 
etry of the mesoscopic object. Here we show that by 
tailoring the geometry of the mesoscopic magnet one can 
produce deterministic dynamic switching between well 
defined degenerate topological states using only in-plane 
magnetic pulses. We present a simple model that ex- 
plains the mechanism of the vortex nucleation and origin 
of robustness of the vortex polarization. Confirmation 
of the model is accomplished by using numerical simula- 
tions. 

Controlled manipulation of magnetic domains in ferro- 
magnet nanostructures have recently opened opportuni- 
ties for novel fast, high-density, and low-power memories 
with novel architectures^^. Any perspective magnetic 
memory architecture, such as sub- micron non- volatile 
magnetic memory, has to have (1) a well-defined switch- 
ing field used to set the memory bits, and (2) a repro- 
ducible switching behavior using simple sequence of ex- 
ternal magnetic field pulses. Therefore, the dynamics of 
the switching between different ground states has to be 
understood in detail. 

Recent advances in fabrication technology at nano- 
scale have enabled studies of magnetic systems that are 
well defined in all three dimensions on a nanometer 
length scale (< 100 nm)^^^. The size reduction of nano- 
magnets leads to novel spin topological states such as 
vortex state, C- state, S-state, flower state, etc.^, and 
to simplified transitions between these states in exter- 
nal magnetic field. The transition between ground states 
in such nanomagnets is of fundamental importance^. It 



is governed by competition between the magneto-static 
energy and exchange energy, and it is influenced by the 
magnetic material used and by the choice of the nano- 
magnet shape. For example, the magnetization of disks 
in zero field can be oriented in-plane, out-of-plane, or 
vortex state can be created, depending on the disk diam- 
eter and thickness^'^. In the disk the flux-closure mag- 
netic state reduces the long range stray fields, i.e. reduces 
the magnetostatic interaction between neighboring disks. 
Therefore, such disk magnetic systems have a potential 
for high-density magnetic storage elements, with bits rep- 
resented by chirality and polarity of a basic vortex state. 

In sub-micron sized disks four possible states have 
to be controlled. It was shown experimentally that in 
the nanomagnets with broken rotational symmetry, chi- 
rality can be controlled easily by in-plane field of se- 
lected direction^^^. At the same time the polarity of 
the vortex core, which represents the second bit, can 
be controlled by an out-of-plane magnetic field^^, spin- 
polarized current ^^ , high-frequency in-plane magnetic 
field^^, or by an in-plane magnetic pulse of precisely de- 
fined amplitude and duration^^. 

In our previous work we have proposed a new prospec- 
tive shape of a nanomagnet with broken symmetry which 
permits control of chirality and polarity bits by the ap- 
plication of the in-plane field only^^. In this paper we 
analyze the mechanisms that establish specific chirality 
and polarity values in such Pacman-like (PL) nanomag- 
net by taking a closer look at energies and dynamics that 
govern the switching processes. We show that the polar- 
ity and the chirality of the vortex core nucleated in the 
decreasing in-plane magnetic field is implicitly defined by 
the direction of the magnetic field with respect to missing 
sector of the PL nanomagnet. 

We consider a magnetic dot of cylindrical shape. To 
construct a PL structure it is necessary to remove an 
outer sector that is 45° wide, and has 1/3 of the disk 
radius (see Fig. [1]). Due to symmetry analysis that will 
be presented later, we choose the orientation of x^ z axes 
such, that they define mirror symmetry plane ay, which 
leaves PL object invariant. Another symmetry operation 
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FIG. 1. Geometry of Pacman-like nanomagnet. The struc- 
ture is symmetric with respect to reflection plane ay. 




is mirroring through the plane x, y noted in following text 

First, we define polarity 7f[/] and chirality x[f] vectors 
as functionals of an arbitrary vector field /: 



X[f]= /'rx(/(f)-7f)dfi 



(1) 



Polarity is just simple average value of the field, while 
chirality resembles the definition of momentum of quan- 
tity / in classical mechanics. Subtraction of polarity in 
the expression for chirality is necessary due to chirality 
invariance with respect to the origin coordinate system 
choice. Mostly we are interested in z component of po- 
larity and chirality. Integration domain Q is over the 
volume of the PL nanomagnet. 

Let's consider that the nanomagnet is placed in a 
strong in-plane magnetic field that has an angle (p with 
the X-axis. To emulate magnetic response by the missing 
sector, we can consider the PL nanomagnet as a super- 
position of a full disk and a set of microscopic magnetic 
moments in a removed sector. These additional moments 
have to have the same value and to be oriented in oppo- 
site direction to the magnetic moments in the disk (Fig. 
[2]). In the first approximation all microscopic moments 
are parallel. Neglecting higher than dipolar moments, 
missing sector behaves as a dipole with a moment mcut 
positioned in the center of mass of the sector fr: 



f Mdn' ^T = 7^ ffdn'. 



(2) 



Minus sign reflects that dipoles of opposite orientation 
have to be added in order to eliminate dipoles in the 
missing sector (see Fig. [2j). The integration is over the 
volume of the missing sector Q^ 

To calculate the magnetic field of PL nanomagnet con- 
sider first the full disk. In magnetic fields exceeding the 
saturation field the magnetic polarization M is parallel 
to the direction of the applied field i?ext throughout the 
disk volume. Internal magnetic field H is also uniformly 
oriented in parallel with M. As small piece of mate- 
rial is removed, all the fields change slightly. To correct 



FIG. 2. (a) The magnetization of the Pacman-like nanomag- 
net is a superposition of the uniform magnetization of a full 
disk (red arrows) and the magnetization of the missing sector 
that is equal and opposite to the one of the full disk, (b) 
The sum of the compensation moments in the sector creates 
a dipole rricut which asymmetrically interacts with the local 
magnetization in the nanomagnet. 



the internal magnetic field, the field of magnetic moment 
^cut given by equation ([2j) has to be added to the orig- 
inally homogeneous internal field of the full disk. The 
removed part thus creates a dipole which induces mag- 
netic field with non-zero chirality x[^]^ if (f 7^ 0°, 180°. 
Since dipoles partially follow the field orientation, the 
non-zero chirality of M is expected. 

Next, we explain qualitatively the mechanism which 
determines vortex core polarity. Taking into considera- 
tion magneto-statics only, the state with positive polarity 
7r[M] is energetically equivalent to the state with nega- 
tive polarity — 7r[M], if no external field in z direction 
is applied. That can be easily seen by writing the total 
energy functional: 




3M- {r-f')]^' ■ (f-f) 



■ Mdn 



dVtdQ.' 



(3) 



Changing the sign of Mz does not change first two in- 
tegrals of the equation (|3]), since these two parts of the 
energy functional are quadratic in Mz and its derivative, 
respectively. The only linear part in Mz is in the third 
integral. But a change of sign of Mz would not influ- 
ence the total energy, because i^ext has no z component. 
Integration domain Vt is not changed by reversing z^ be- 
cause symmetry operation a^ - reflection from the plane 
xy transforms PL nanomagnet to itself. 

The final state of vortex core polarization is deter- 
mined by magnetization dynamics. Time evolution of 
magnetization is described by phenomenological Landau- 
Lifshitz- Gilbert (LLG) equation: 



dM 



-jM X Hef^ ^a[ M X 



dM 



(4) 




FIG. 3. Angular dependence of the vortex nucleation pro- 
cess. Top: angular dependence of vortex nucleation field. Di- 
rection of the nucleated vortex polarity and chirality are also 
indicated for each quadrant. Bottom: two different processes 
of vortex nucleation depending on initial magnetization direc- 
tion with respect to PL's symmetry plane {^ — 15 o and ^p — 
75o): from uniform magnetization state the magnetization 
transitions to S-shape (dots) or C-shape (stars) configura- 
tions and equilibrates to a vortex state with specific polarity 
and chirality. 



Here we show that this equation itself contains polarity 
symmetry breaking mechanism. First, we apply strong 
external in-plane field in direction that has an angle Lp 
with the X-axis (Figl2]). Then we slowly (adiabatically) 
decrease the external field amplitude to the level that is 
just above the vortex nucleation field. By adiabatic field 
change we mean that the change of the external field with 
time is so slow, that the energy dissipation keeps the sys- 
tem very close to the local minimum at all times. Then 
^^ ~ everywhere. Any dynamics means also dissipa- 
tion of energy due to term that is proportional to a in 
equation (|4j). Therefore, at local minima we also have 
^eff ^ 0, otherwise it is not possible to satisfy equa- 
tion (HI) with vanishing left side. Now, having i?ext just 
above nucleation field, we decrease the external field by 
a small value A J?. To first approximation the effective 
field is Hq^ = AH. When looking at the dynamics of 
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FIG. 4. Angular dependence of the in-plane vortex nucle- 
ation field Bnuc and threshold out-of-plane field Bz necessary 
to reverse the polarity of the entering vortex. 



local magnetic moments shortly after the external field 
is decreased, we can neglect the damping term in the 
LLG equation, since a <C 1. The dynamics of Mz is then 
governed by the equation 



dt 



= -jiM^AHy-MyAH^). 



(5) 



The right-hand side of the equation (J5j) is not zero lo- 
cally, nor it is in average, due to asymmetry of the PL 
nanomagnet with respect to the direction of the applied 
field (if (/? 7^ 0°, 180°). If the actual value of the external 
field is lower than the nucleation field, the non-zero po- 
larity of the nanomagnet will emerge. Paths towards two 
equivalent minima characterized by iM;^ are energeti- 
cally equivalent, but change of external field with time 
is biasing the time evolution of magnetization towards 
configuration given by eq. ([5]). 

Confirmation of the above model approximation can be 
obtained by numerical simulations. We have performed 
numerical simulations of PL nanomagnet using OOMMF 
software package^^. Parameters used in the simulation 
are: outer radius R = 3b nm, thickness (in z-direction 
- not shown) h = 40nm. Material used in the calcula- 
tions was Permalloy Ni8oFe2o (Py), with following mate- 
rial parameters: exchange constant A = 13 x 10~^^J/m 
and saturated magnetization Mg = 8.6 x lO^A/m. 

In Figure [3] we show the dependence of the nucleation- 
field amplitude on the applied field direction. The nu- 
cleation field is defined as the applied magnetic field at 
which the non-zero z component of the magnetization 
polarity 7r(M) appears. External field is adiabatically 
decreasing from 150 mT to zero in selected direction. By 
adiabatic change we mean repeated process of decreasing 
field by 2 mT step, followed by full system relaxation. 

We would like to note the symmetry properties of PL 
nanomagnet. In Fig. |3l(top) each quadrant corresponds 



to a vortex ground state of the nanomagnet with specific 
chirality and polarity shown in corresponding corners. 
The nanomagnet relaxes into that remanent state from a 
uniform magnetization along an angle within the specific 
quadrant. As can be seen in FiglH angular dependence of 
the nucleation field can be reconstructed from the depen- 
dence for (/p G [0; 1^] by inversion and refiection through 
xz plane ay . Inversion symmetry of the graph shown in 
Fig|3] is the consequence of the time-reversal symmetry. 
The refiection symmetry with respect to xz plane shown 
in Fig|3]is related to refiection from (jy - geometrical op- 
eration that transforms PL nanomagnet to itself. 

The simulation results show the existence of two dis- 
tinct vortex core nucleation regimes (Figl3] (bottom)). 
For large angles {cp G [50°; 90°]) - vortex nucleates from 
C-state magnetization pattern. This form of nucleation 
is not robust in the sense that even small out-of-plane 
field Bz — ImT is sufficient to alter the resulting polar- 
ity of the nucleated vortex along the direction of applied 
field Bz. Instead, for smah angles {cp e [0°;48°]) - the 
vortex nucleation path is different. Just above the vor- 
tex nucleation field, the magnetization of PL nanomagnet 
forms an S-state. This configuration consists of two re- 
gions with opposite signs of the curvature of field lines^^. 
Meanwhile, there is only one curvature of field lines just 
below the vortex nucleation field. The process of vortex 
core nucleation in this case involves reversal of magnetic 
moments in a part of the nanomagnet. This reversal pro- 
cess proceeds through an out-of-plane motion of the local 
magnetic moments, resulting in robust vortex core polar- 
ization despite the presence of small external fields in z 
direction. In figure |4] we show the angular dependence of 
the maximum external field B^ for which the PL nano- 
magnet is able to sustain nucleation of vortex polarity 
opposite to the direction of applied external field. 

The C and S-shapes of magnetization can be explained 
by the position of perturbing dipole mcut- As can be 
seen from the Fig|2l PL nanomagnet is divided into two 
domains with opposite sign of magnetic field circulation 
generated by mcut- The size of this two domains are 



determined by orientation of mcut- However, exchange 
interaction tends to align local moments in parallel, thus 
there exists a critical angle (around 48° in our geome- 
try) , beyond which the region with minor curvature does 
not exist. Detailed energy balance between exchange and 
cavity (sector) demagnetization determines the scenario 
of vortex nucleation and its eventual robustness with re- 
spect to external perturbation. 

To summarize, in this work we provide simple argu- 
ments that elucidate the origin of driving mechanisms 
for nucleation of magnetic vortex with controlled chiral- 
ity and polarity. We have also found the regime of PL 
nanomagnet operation in which the final vortex state is 
independent on weak disturbing external field. This is a 
promising finding to consider if using PL nanomagnet as 
a memory element in bit-patterned media or as a gener- 
ator of magnetic vortices of desired polarity and chiral- 
ity for microwave applications. Weak interaction among 
the elements as well as robustness to small external field 
perturbations makes PL nanomagnet very suitable for 
operation. 

Finally we note, that the sub- 100 nm PL nanomagnet 
is not the only unique design offering control of chiral- 
ity and polarity by in-plane magnetic field. Qualitatively 
similar results are obtained in simulations for different 
sizes and shapes of the missing sector. According to our 
simple model, the necessary ingredients are symmetry of 
the object, the demagnetization field strength of removed 
part and shape anisotropy induced by removed part. Ro- 
bustness of vortex polarity against Bz is based on vortex 
- anti- vortex annihilation during vortex core nucleation. 
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